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Abstract 



Q_i| Here we study the effect of the non-minimal coupling j^€ flua d l 'A a on the static 

potential in multiflavor QED3. Both cases of four and two components fermions 
are studied separately at leading order in the 1/N expansion. Although a non- 
local Chern-Simons term appears, in the four components case the photon is still 
massless leading to a confining logarithmic potential similar to the classical one. 
In the two components cs.se, cis expected, the parity breaking fermion mass term 
generates a traditional Chern-Simons term which makes the photon massive and we 
have a screening potential which vanishes at large inter-charge distance. The extra 
non-minimal couplings have no important influence on the static potential at large 
inter-charge distances. However, interesting effects show up at finite distances. In 
particular, for strong enough non-minimal coupling we may have a new massive pole 
in the photon propagator while in the opposite limit there may be no poles at all 
in the irreducible case. We also found that, in general, the non-minimal couplings 
lead to a finite range repulsive force between charges of opposite signs. 
PACS-No.: 11.15.Bt , 11.15.-q 
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1 Introduction 



One of the most intriguing and longstanding problems in high energy physics is a com- 
plete understanding of the mechanism of color confinement in 4D QCD. Several models 
and techniques have been used to gain a deeper insight on this problem. In particu- 
lar, supersymmetry has been used in the four dimensional theory treated in 1] while 
bosonization has played an important role in the analysis of screening/confinement is- 
sues in D = 2 models El El- Those models have in common that it is necessary to 
have massive fermions to produce confinement. This happens also in parity preserving 
QED in 2 + 1 dimensions [5] which will be treated in this work altogether with the parity 
breaking (2 components fermions) case. Technically, although bosonization is not so well 
developed in D > 2 , we can always integrate perturbatively over the fermionic fields (see 
jo] I7[ IH1 1§1 ITm ITTj land derive an effective bosonic action for the vector field. Such effective 
action is in general non-local even in D = 2 From this bosonic action we have an 
expression for the vector boson propagator including vacuum polarization or even higher 
order effects. A detailed analysis of the analytic structure of the propagator already re- 
veals, qualitatively, the large distance behavior of the interaction. In fact, even in other 
non-equivalent approaches to study confinement, like Schwinger- Dyson equations, see e.g. 

the infrared properties of the vector field propagator including vacuum polarization 
corrections play a key role. Here we formally minimize the effective action and calculate 
the energy V(L) between two static charges separated by a distance L by solving the 
resulting differential equation. Basically, this is the route followed in ^21 E31 E] in the 
case of QED3 with two components fermions (irreducible QED). In that case an usual, 
linear in momentum, Chern-Simons term is generated which makes the photon massive 
[To] leading to a screening potential in opposition to what happens classically. In ^H] we 
have added a quartic interaction of Thirring type Lth = g 2 ^^^ 1 ) 2 /N and checked that 
the screening scenario prevails again in multiflavor irreducible QED 3 , at least at leading 
order in 1/N, see also |T7J [IB]. By passing, we noticed in JE] that in the case of four 
components fermions with a parity symmetric mass term, henceforth called reducible 
QED, the vacuum polarization is not strong enough to change the classical picture and 
the potential remains logarithmic confining at leading order in 1/N. In all the works 
[T2*[ IT3[ ITU ITT)] the fermions were minimally coupled to the eletromagnetic potential. The 
aim of this work is to analyze the effect on the screening/confining scenario of adding a 
non-minimal coupling term which preserves gauge symmetry and it is rather natural in 
2 + 1 dimensions, namely, we add F^j^ = e^ va d v 'A a j fl to the minimally coupled QED 
action. In the irreducible case this extra coupling corresponds to a magnetic moment 
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interaction of the Pauli type. Such non-minimal coupling has been considere before in 
[El 1201 I2U 1221 123 121] in different contexts. It seems to lead to anyons without the need 
of a topological Chern-Simons term as argued in |2~T1 121] and [22] ( see however |2H])- 
Our interest lies in the fact that F^j^ breaks parity explicitly and the presence of parity 
breaking terms is very important for the pole structure of the photon propagator. 
We start with the partition function: 
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T>A fl Y\ T>ip r T>ip r exp 
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Summation over repeated flavor index r (r = 1, 2,...,N) is as- 



sumed. The constant /i represents the magnetic moment of the static charge while y sets 
the strength of the dynamical fermions non-minimal coupling. Those couplings have mass 
dimension —1 and —1/2 respectively. The external current is generated by a static charge 
at the point (xi,x 2 ) = (L/2,0), i.e., 



jl xt = QS(x 2 )S(x 1 - -)S 



I/O 



(2) 



Integrating over the fermions we have an effective action for the photons: 



J VA^exp \yi j d 3 x 
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eA^{k)+yF^(k) 
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The logarithm has been evaluated perturbatively in 1/JV. Due to Furry's theorem, only 
even number of vertices contribute. Since each vertex is of order 1/yN the leading 
contribution with two vertices, that corresponds to the vacuum polarization diagram, 
will be iV-independent due to the trace over the internal femion lines. The next to 
leading contribution with four vertices is of order 1/N and will be neglected henceforth. 
In particular, our results will be exact for N — > 00. The quantities A^(k) and F^(k) 
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represent the Fourier transformations of A^(x) and F^(x) respectively and IP" is the 
polarization tensor: 



d 3 p 

J2^y 



tr 



1 



-r 



-7 



(4) 



ft — m + ie (ft + If) — m + ie 
The calculation of IP 1 ' is regularization dependent. Using dimensional regularization 
which preserves gauge symmetry and does not add any artificial parity breaking term, we 
have : 



n^u/(/c) — m kfjzj) + -^e^ VQ k 



(5) 



Defining 2 = A; 2 / (4m 2 ), in the case of two components fermions, in the range < z < 1, 
we have 



/i 
/a 



1 . 1 + 
— f= In f= 

^ 1- V* 

i + ^A 

^ \ 2 17 



(6) 



While for z < 



/ 2 



arctan v — 2; 



If 1- 



-/1 



(7) 



Above the pair creation threshold z > 1 the effective action will develop an imaginary 
part which we are not interested in and do not write it down here. In the simpler case of 
four components fermions we have instead: fx — > , f% — > If 2 ( fi — * , J2 — > 2/ 2 )• 
From (JHJ) we can read off the effective action: 



5, 



cff 



A^(x) ^4^ A y (y) + (A v + nF v )3U®{x - y) 



Where D*(x,y) is the inverse of the propagator in the space-time. Minimizing the 
effective action we deduce: 



Mv) 



d 3 xD Pa (y,x) [& t (x) + Ai£ y 0"j£ t (s)] 
d 3 k 



D Pa {k) / d»xe«-(»-) [ 3 : xt (x) +f ie^jl t (x)} 



(9) 
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Introducing the dimensionless couplings: 



Ci = e 2 / 16iim 
c 2 = e^/An 

We can write the propagator in momentum space as: 



D^(k) = ag^ + b yg^ 
For < z < 1 the coefficients are given by 



£2 I ^ 



(10) 

(11) 
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a + b 

c 
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while for z < 0, 



4m 2 (z 
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8m 2 ^D + D_ 
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16m 3 ;? D + D- 
(ci ± \fzc 2 ) 2 G± - v^ci 



G± = v^/ 2 ± /i 
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ci(c 2 + zc 2 .) 2 



4 2 + zB 2 
^c 2 - ^)/i - 2c x c 2 f 2 z 



lQm 3 z(A 2 + zB 2 ) 
(c 2 + zcj) 2 f 1 -2zcjc 2 ; B 



ci(4 - zc l) 



(17) 



(19) 
(20) 



Returning to the calculation of ©, since the external current is time independent the 
integral J dxoe~ tk ° Xo = 2^5{ko) allows us to exactly integrate over ko which amounts 
to set k^k^ = —k\ — k\ = —k 2 < inside ©• We can integrate over dx\dx 2 using 
5(x\ — L/2)5(x 2 ). The integral over the angle part of d 2 k = kdkdO gives rise to a Jo Bessel 
function. Thus, we are left with the integral over k = \Jk\ + k\. Placing the negative 
charge — Q at (xi, x 2 ) = (— L/2,0) we finally have the energy of the pair separated by a 
distance L: 
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V(L) = -QA ( yi = -L/2,y 2 = 0) 

= J dkk {a + b + ifidk 2 ^ J (kL) (21) 

The last formula shows that the non-minimal /x-term can only contribute to V(L) if the 
photon propagator contains a parity breaking piece (c ^ 0). 

In the next two sections we split the discussion into the cases of four (reducible) and 
two (irreducible) components fermions. 

2 Reducible QED (4x4 representation ) 
2.1 Effective action and pole analysis 

In this case, the inclusion of vacuum polarization effects leads to the following non-local 
Lagrangian density: 



C 



/IV 
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+ (A v + ^F v ) jl 



cxt 



:22i 



Where f 2 = f2{—0/4m 2 ) is given in © and (JJ|). Notice that, within dimensional regular- 
ization IL^ is parity symmetric (/i = ) but still, as a consequence of the parity breaking 
non-minimal coupling there appears a non-local Chern-Simons term in ()22j) which breaks 
parity . 

Concerning the pole structure of the propagator, setting f 1 — > and f 2 — > 2f 2 in (fTBj) 
we obtain D + D_ = z [2(ci + \fzc 2 ) 2 f 2 — C\] [2(ci — \^zc 2 ) 2 f 2 — Ci]. Since f 2 < —4/3 and 
c\ > it is clear that the only pole we have appears at the origin z = . Therefore, quite 
surprisingly, the non-local Chern-Simons term in f!22l) is not able to make the photon 
massive like its local counterpart. Although parity is broken, the photon is still massless. 



2.2 V(L) 

For the calculation of V(L) it is important to find the singularities of (j!8j) and |TT)|) . By 
inspection we observe that the denominators of a + b and c can only vanish either for z = 
or A = = B. However, since fi = it is clear from (|20|) that the later possibility also 
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requires z = 0. Thus, we conclude that the expression to be integrated in V(L) contains 
only one singular point which is a simple pole at z = 0. By passing, this implies that we 
have no tachyons in the reducible case. 

Substituting (JTHjl and (fl^jl. with /i — ► and j 2 — > 2/ 2 , in (EH) we derive, after trivial 
manipulations: 

r°° dk 

V(L) = J q T J (kL)G(k) (23) 
Where, besides the Bessel function J , we have 

C(k) = ® 2 ° l Cl<yl - 2Cl ^ + 2c2Z ^° 2 - 2/iCl ^ (24) 
" " 2vr (2^/ 2 + c 1 ) 2 -4c?/ 2 (l-c 1 / 2 ) 

Note that / 2 < —4 / 3 which confirms that k = is the only singularity in (|2*3*j) . The way it 
stands, the integral in does not exist since Jo(0) = 1 and G(0) = — 3Q 2 /[27r(3 + 8ci)] 
are both finite and non- vanishing . In order to eliminate the infrared divergence at k — 
we make a subtraction, namely : 

r°° (iu 

V(L) - V(L ) = lim / — [J (kL) - J (kL )}G(k) (25) 

We can recover the classical result by taking c\ = = c 2 which gives G(/c) = —Q 2 /27r. 
In this case the integral ()25j) can be easily calculated furnishing the well known confining 
potential: 

V(L)-V(L ) = ^lnA (26) 

However, in the general case we have to calculate the integral numerically One exception 
is the limit m — > oo. First, notice that if we simply take m — > oo the effective action 
122(1 becomes the classical Maxwell action at leading order and all quantum information 
is lost. On the other hand, if we take m — > oo while keeping c\ and c 2 finite the effective 
action (1221) becomes: 



£e// = ~\ (l + ^) V" + + + ^) iext + ^ (^) (27) 

Correspondingly, this amounts to substitute G(k) by G(0) in ([25)1 . Therefore we have : 

30 2 

^(L)-V(Lo) = -G(0)ln(L/L ) = p^^log(L/L ) 

= ^log(L/L ) (28) 
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Figure 1: The static potential for pure QED 3 (reducible) with (ci,C2,/x) = (1,0,0) and 
m — 1 (dark dots), m — 0.1 (light dots). The dark solid curve corresponds to the m — > oo 
result (3/11) InL and the light solid one to the fit (3/11) In L - 0.48/L + 0.046/L 2 + 0.44. 



Where 

Qscr. = ^1/2 ( 29 ) 

L 6-7rmJ 

Thus, we see that at leading order in 1/m the vacuum polarization has a mild effect on the 
classical static potential. It leads to a screening of the static charges but it is not strong 
enough to change the confining nature of the potential. In particular, the non-minimal 
couplings introduced by the constants 7 and /i have no influence at all. In fact, this is 
not surprising since we are keeping c\ and c 2 fixed and taking m — > 00 which corresponds 
to e ~ \Jm and 7 ~ 1/y/m. This is a situation where the minimal coupling certainly 
prevails against the non-minimal one. 

In order to clearly outline the effect of the non-minimal coupling we must keep all 
parameters (ci, C2, m) finite and calculate V(L) numerically. Our results for the reducible 
case are shown 1 in figure 1 (pure QED) and figure 2 (non-minimally coupled QED). Fixing 
q = 1 we see that in both figures the value m — 1 is already very close to the m — > 00 
case. This could also be interpreted as a confirmation of our numerical integration which 
is rather trick due to the oscillations of the Bessel function. The remarkable effect of the 
non- minimal coupling is a repulsive force in a finite range < L < L m i n . Numerically we 
have noticed that L min ~ 1/m. In particular, in figure 2 L min ~ 2.5 and 0.25 respectively 
for m — 0.1 and m — 1. The effect produced by c 2 7^ and fi = can be similarly caused 

1 In figures 1 and 2 the symbol V stands for the difference V(L) — V(Lq — I) while in figures 3-7 it 
represents V(L). The potential is always depicted in units of Q 2 /2tt. 
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Figure 2: The static potential for QED3 (reducible) with non-minimal coupling 
(ci,C2,/i) = (1,0.5,0) and m — 1 (dark dots), m = 0.1 (light dots). The dark solid 
curve corresponds to the m — > 00 result (3/ll)Log[L] and the light solid one to the fit 
(3/11) InL - 0.42. 



by making /i 7^ and C2 = 0. Furthermore, they can also compete such that the repulsive 
force can be turned off if we set \i and C2 accordingly. The same happens in the irreducible 
case (next section) and the figures we have are similar (for small L) to the figures 3-5. 
Therefore we skip them. The logarithmic fittings in figures 1 and 2 are in agreement 
with the analytic result derived in the appendix of [SHI for pure QED (reducible) , namely, 
V(L) = —G(0) InL + constant + h(L) where h(L) should fall off at least as fast as 1/L 
as L — > 00. In summary, the new couplings 7 and \x can only play a role at finite L. The 
point is that if we use a 4 x 4 representation only a higher order Chern-Simons term will 
be generated due to the parity breaking non-minimal coupling but as k — > this term 
will be negligible if compared to a Maxwell term and the photon remains massless as it is 
classically. Consequently, the origin k = will dominate the calculation of V(L). Around 
that region the higher momenta coupling terms which multiply 7 and fi can be dropped. 

3 Irreducible QED (2x2 representation ) 
3.1 Effective action and pole analysis 

Due to the parity breaking term of the effective action now is more complicated: 



eff 



e7 D/ 2 + f e 2 - 4 7 2 D ] ^ 



87rm 



16tt 



d u A c 
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[IV 



-4 7 2 D\ e 7 ' 
~T7. /2 — — n 

lonm I 



2tt' 



+ (A u + iuF u ) & 



ext. 



(30) 



Where f 2 = / 2 (— D/4m 2 ) and fi = D/4m 2 ) are given in (JHJ) and ((7J) . 

For finite fermion mass the analysis of the singularities of the photon propagator is 
much more involved. First of all, the pole at z = in (fl3jl .which also appears in the local 
Maxwell Chern-Simons theory, is not a physical one. In particular, it disappears if we 
look at the gauge invariant propagator < F )1 (k)F u (—k) >. Thus, the physical poles can 
only come from either D + = or D- = 0. Since G + < —2 the first possibility is ruled out, 
see (Jinj). Therefore, we have to examine the equation _D„ = (ci — A/ic 2 ) 2 G_ — \fzc\ = 
with G- being a monotonically decreasing function in the range 1 < G_ < 2. We have 
not been able to find an analytical solution to that equation but we have made a rather 
detailed analysis on the number of poles (n) in the region < k 2 < Am 2 according to 
different coupling values. One can show analytically that, 



c 2 > ci + ^fc[ — > n = 2 (31) 

c 2 < -1/4 -> n = (32) 

-1/4 < c 2 < ci - Jc[ -> ti = fl if ci>l (33) 

ci — \fo[ < c 2 < c\ + a/ci — > n = 1 if ci > 1 (34) 

< c 2 < ci + ^/cT — > n = 1 if ci < 1 (35) 



The missing region ci < 1 and —1/4 < c 2 < is rather awkward and only numerical 
results have been obtained. In particular, besides n = 0, 1, 2 poles we also found it possible 
to have n = 3 poles if the QED coupling is very small < c\ < 0.041. 

It is remarkable that for c 2 > 0, the poles follow basically the same pattern ^H] of 
QED with a Thirring term Crh — g 2 (ip t 'j fJ "ip 1 ) 2 /N ^H] with the Pauli coupling 7 playing 
the role of g. Namely, for strong enough Pauli coupling, when compared to the QED 
coupling , a new pole appears besides the one we have in pure QED. If 7 is reduced this 
extra pole disappears. Besides, for very strong QED coupling we may have no real poles 
at all. That also happens in pure QED as we remarked in JH]- It is tempting to ascribe 
those similarities with the QED plus Thirring case to the fact that both couplings g and 
7 have the same mass dimension mT 1 ! 2 . 

Finally, we stress that it is impossible to have a massless photon in the 2x2 repre- 
sentation whatever value we choose for c\ and c 2 . 
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3.2 V(L) 

Now the integral we need to evaluate is: 



poo 

V(L) - V(L Q ) = lim / dkk[J (kL) - J (kL )]H(k) 

Jx 



(36) 



Where, 



H(k) 



Q 2 a (c\ + zc\) 2 [(c{ - zc\ + 2jimc 1 c 2 z)f 2 + 2c x c 2 fi - //m(cj - zc\)fi - a] 



8vrm 2 [( c 2 + zc ^2~ h _ 2zc j C2 ]2 + ^[( c 2 + zc 2 )2 / 2 _ Cl(c 2 _ zc 2 )]2 

(37) 

In the reducible case the integral for V(L) was dominated by the pole at the origin 
k = \/kf + k\ = 0. Now, in order to have a singularity both expressions inside brackets in 
the denominator of (}3Tj) must vanish at the same time. This will never happen at the origin 
since /i(0) = —2 and f 2 (0) = —4/3. A detailed analysis reveals that this can only happen 
for some real k > if we fine tune c\ and c 2 . The fine tuning is only possible in the small 
region —1/4 < c 2 < and < c\ < 0.1354. Since k^k^ = —k\ — k\ < the singularity is 
interpreted as a tachyonic pole. Henceforth, we exclude the above fine tunning from the 
parameters space which guarantees that we are tachyon free. Consequently, there will be 
no singularity in the integration range of (|36|) and we need to calculate it numerically. 
An important consequence of the absence of singularities is that we can interchange the 
limit and the integral: lim^oo J °° dk k Jo(kL) H{k) = J °° dk lim^oo k J Q {kL) H{k) = 0. 
Therefore, we certainly have a screening potential in this irreducible case no matter what 
we choose for the couple (ci, c 2 ) or for the external charge magnetic moment /i. This result 
has been confirmed by our numerical calculations using MATHEMATICA software. For 
all cases the potential tends to zero as L — ► oo. The specific details of the potential are 
exhibited in the figures. 

First, we show in figure 3 the effect of the external charges magnetic moment fj, in the 
pure QED case without the Pauli term (c 2 = 0). It is interesting to notice that for large 
values of \x a repulsive force appears for finite L which changes the form of the potential 
compared to what one has in previous works of pure QED ^21 E3 E] • 

In figure 4 we take /i = and analyze the effect of the magnetic moment of dynamical 
fermions ( Pauli term ) at fixed QED coupling C\ — 1. Similarly to the reducible case, the 
Pauli term leads also to a new, if compared to pure QED, repulsive force which is placed 
at small values of L. That happens even for tiny values of the Pauli coupling c 2 . We 
have also let the couplings /i and c 2 appear altogether, the final output is shown in figure 
5. It is possible to turn the repulsive force which appears for a fixed value of c 2 into an 
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Figure 3: The static potential V(L) for QED 3 (irreducible) with different values of the 
static charges magnetic moment \x = —2, 0, 2, 20 respectively from the darkest to the 
lightest curve. All curves were obtained for m — 0.1 , c± — 1 and vanishing non-minimal 
coupling C2 = 0. 
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Figure 4: The static potential V(L) for QED 3 (irreducible) with fj, — and different 
non-minimal couplings c2 = 0.5, 0.2, —0.5 respectively from the darkest to the lightest 
curve. All curves were obtained for m — 0.1 , cl = 1. 
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Figure 5: The static potential V(L) for QED 3 (irreducible) with /i = 3, 0, —4 respectively 
from the darkest to the lightest curve and (ci,c 2 ,m) = (1,0.5,0.1). 



attractive one by changing the static charges parameter [i. Due to a term that depends 
on the product /ic 2 , see (J27J), now the influence of /i is opposite to the pure QED (c 2 = 0) 
case. As we increase [i the repulsive force diminishes. In the specific case of fig.3, i.e., 
(ci,C2,m) = (1,0.5,0.1) we checked that no repulsive force appears for \i > 1.67. 

3.3 Large mass limit 

If we repeat the procedure of last section and send m —>■ oo while keeping c\ and c 2 
constant, the effective action (J3TIJ) blows up at leading order. So we found more useful 
instead to take m — > oo and keep e and 7 fixed in order to have a local theory, which 
becomes : 

a , C(<V^) 2 



Ceff = -(l + ^)^F^-i- e ^(e 2 -4 7 2 D)^ + 

+ (A u + ^F u )j^ t + o(-) (3* 



m r 

Besides the well known generation of a Chern-Simons term, we also have a charge renor- 
malization : 

Qren — ~ 7772 ' (39) 

(1 + f) ' 

which is due to a cooperative effect with the Pauli 2 interaction which in its turn also 
generates a higher order Chern-Simons term. This higher derivative term leads in general 
2 The 2x2 Dirac matrices satisfy the algebra [7^,7^] = —2ie^ ua j a . Consequently, the non-minimal 
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to an extra pole in the photon propagator. Both poles are real for c2 > —1/8 (e^/n > 
—1/2). The formula (J21j) becomes in this case: 

v <n _o^f^±nY\ r dkkMkL) 

m ^°° l 1 ~ 2tt V 7 4 J Jo (k 2 + m 2 + )(k 2 + m 2 _) 1 } 

Where 



We can check that as 7 — > , which is the case of pure QED 3 , 7 4 (fc 2 + m 2 + ) — > n 2 
while m_ — > e 2 /47r. Thus, we are left with the well known massive photon of the Maxwell- 
Chern-Simons theory which is the effective action for pure QED 3 (irreducible) at m — > 00. 
The simplicity of the photon propagator allows us to calculate V(L) analytically. We just 
need the identity l/{k 2 + m 2 _) — l/(k 2 + m 2 + ) = (m 2 + — m 2 _)/[(k 2 + m 2 + ){k 2 + m 2 )]. 
Consequently, 

V m -*ao(L) = %- * [K (m + L) - iT (m_L)] (42) 

27T / 1 + 2ej 



For 7 — > we reproduce, at leading order in 1/m, the result of pure QED 3 obtained in 
|12j . Due to the relative sign among the modified Bessel functions the potential due to 
two different poles is now bounded for L — > 0, since Kq(x — » 0) — > — ln(x), and (|4^j) 
becomes constant as L — > (figure 7). Thus, the attractive force vanishes as L — » 0. 
For finite fermion mass the vanishing force can turn into a repulsive one. In figure 6 we 
compare the potential Vm^oo with the numerical results obtained for m = 0.1, 1 in units 
such that e 2 /(167r) = 0.1 and t^jn = 0.5. 

The agreement, for most of the values of L, of the analytic formula ()42|) obtained at 
leading 1/m approximation with the numerical results already at m = 1 is impressive. Of 
course, we have checked that larger values m > 1 lead to an even better agreement. The 
repulsive force close to the origin only appears in the presence of a non-minimal coupling 
7 7^ ( or fi 7^ ) and a finite fermion mass. 

coupling term can be interpreted in the irreducible case as a magnetic moment interaction of Pauli type: 
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ure 6: The static potential V(L) for QED 3 (irreducible) for m — > oo (solid line), 
= 1 (dark dotts) and m = 0.1 (light dotts) with couplings mc\ = e 2 /(167r) = 0.1 and 
= e^ijix = 0.5. 
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Figure 7: Blowing up the small L region of figure 6. 



15 



4 Conclusion 



In this work we have calculated, at leading order in 1/N, the influence of a non-minimal 
coupling term on the static potential of QED with four and two components fermions. In 
the four components case the only effect of the vacuum polarization at large inter-charge 
separations is a screening of the original static charges. The potential keeps its classical 
shape V(L) = — (Q^ en /27i)\iaL at L — ► oo. Although parity is broken due to the non- 
minimal coupling the potential is still of confining type. The interesting point is that the 
generated Chern-Simons term is of higher order in the momentum and becomes negligible 
when compared with the Maxwell term as k — > 0. Consequently, the photon does not 
acquire mass which leads to a long range confining potential. Technically, the massless 
pole changes the factor kJo(kL) in the integral for the potential V(L) into Jo{kL)/k. 
Since the Bessel function is an oscillating function with decreasing amplitude, as L — > oo 
the integral is dominated by the region around k — 0. In summary, the extra non-minimal 
couplings play no role in V(L) for large distances as far as our 1/N approximation is valid. 

In the irreducible case of pure QED it is well known that a traditional (linear in 
momentum) Chern-Simons term is generated and the classical pole at k = is replaced 
by a massive pole. This picture remains correct whatever values we choose for the new 
couplings 7 and \i. Therefore we expect a short range screening potential as in pure 
irreducible QED. Indeed, we have obtained numerically V(L — > oo) — > as one can check 
in our figures of last section. Due to the lack of the pole at k — one can take the limit 
L — > oo before the integral in (|25jl is performed and the result will vanish as a consequence 
of Jo(x — ► oo) — ► 0. 

Concerning the presence of massive poles we have found basically three different re- 
gions in the coupling space. For strong enough Pauli coupling, when compared to the 
QED coupling (c 2 > c\ + y/cT), we have a second massive pole besides the one generated 
by the usual Chern-Simons term. Our local limit m — > oo indicates that it is due to a 
higher order Chern-Simons term. If we decrease the Pauli coupling we reach an interme- 
diate region and return to just one pole region as in pure QED. The most embarrassing 
result appears for very weak Pauli coupling ( or very strong QED coupling ), namely, we 
may have no poles at all. That happens, for instance, if c 2 < c\ — ^fc{ and c\ > 1. In 
particular, it occurs in pure multiflavor QED (c 2 = 0) if c\ = e 2 /(167rm) > 1. Recalling 
that our effective action for the photon becomes exact for iV — > oo the no pole region 
might indicate a breakdown of the 1/N expansion in the specific region c 2 < C\ — ^fc{ and 

(a > l). 

Although at large distances the new couplings 7 and \i have not led us to new physical 
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effects, at small distances we found a repulsive force of finite range between charges of 
opposite signs. This effect only appears for a finite fermion mass. As we increase the mass 
of the fermion the potential bents down and the repulsive effect disappears in agreement 
with our effective action (jHHj) at m — > oo. We do not have any explanation for the repulsive 
effect but we intent to return to that question by studying the existence of bound states 
in the future. We should remark that a repulsive effect between charges of opposite signs, 
similar to a centrifugal barrier, was observed before in jTTj in irreducible QED% with a 
Thirring term with positive coefficient (opposite sign to the one used in (THj). However, we 
stress that our numerical calculations were obtained without any approximation for large 
masses or small couplings as in The similarity between the addition of the Thirring 
term and the non-minimal coupling might have its root in the fact that both the Thirring 
and the non-minimal coupling constants 7 used here have the same mass dimensionality 
m -1 / 2 . Studies of duality in vector models in D = 2 + 1 coupled to fermionic matter 
|251 12H1 EZj indicate that there might be a more direct relation between the Pauli and the 
Thirring term in d — 2 + 1 dimensions. 

Finally, we point out that another motivation to study the effect of a non-minimal 
coupling of the Pauli-type is the fact that this term is radioactively generated anyway 
(irreducible case), at least in a Maxwell Chern-Simons theory [SHIES]- Furthermore, it 
has been recently shown in (SHj that if the Chern-Simons term is not present from the 
start the Pauli term will be generated with an infinite coefficient which points to the 
need of having this term from the beginning in order to have an infrared finite theory. In 
our calculations we suppressed a possible Chern-Simons term in the starting Lagrangian 
because we would like to single out the effect of the Pauli-type term. It is clearly desirable 
to have both terms from the start but the profusion of coupling constants make the analysis 
made here much more involved which is out of the scope of this work. 
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